Dynamical creation of entangled bosonic states in a double well 



K. Ziegler 

Institut fiir Physik, Universitdt Augsburg, D-86135 Augsburg, Germany 
(Dated: January 15, 2013) 

We study the creation of a bosonic NOON state from the evolution of a Fock state in a double 
well. While noninteracting bosons disappear quickly in the Hilbert space, the evolution under the 
influence of a Bose-Hubbard Hamiltonian is much more restricted. This restriction is caused by 
the fragmentation of the spectrum into a high-energy part with doubly degenerate levels and a 
nondegenerate low-energy part. This degeneracy suppresses transitions to states of the high-energy 
part of the spectrum. At a moderate interaction strength this effect supports strongly the dynamical 
formation of a NOON state. The NOON state is suppressed in an asymmetric double well, where the 
double degeneracy is absent. 

PACS numbers: 03.65.Aa, 03.65.Fd, 03.67.Bg 

I. INTRODUCTION 

Recent experiments on ultracold gases in optical potentials I1h3(| and experiments on photons in mi- 
crowave cavities [1, Q have demonstrated that it is possible to prepare a Fock state as a pure state 
in a finite-dimensional system. After the preparation of the Fock state, the parameters of the system 
can suddenly be changed (performing a "quench" ) such that the Fock state is not an eigenstate of the 
new system Hamiltonian H . Then the evolution of the many-body state due to the evolution operator 
exp{—iHt) will lead to a random walk inside the available Hilbert space. The visited states include other 
Fock states as well as superpositions of Fock states. Typical questions in this context are: what is the 
probability for visiting different states and how is this affected by the interaction of the particles? A 
natural quantity for measuring this probability is the spectral density function of the Hamiltonian H 
with respect to the initial Fock state [a-Q. 

A classical candidate for modeling the evolution of a Fock state is the Hubbard model [l^ . The 
corresponding discrete Hamiltonian describes the tunneling of a particle between neighboring potential 
wells and a local particle-particle interaction. The Hubbard model for bosons (Bose-Hubbard model) 
was realized as an atomic system in an optical lattice illj. A possible realization of the Bose-Hubbard 
model by photons in coupled microwave cavities was proposed recently by Hartmann et al. An 
anharmonicity of the microwave cavities plays the role of the photon-photon interaction 0] . 

The simplest system for discussing the evolution of a Fock state within the Hubbard model is a double 
well, where particles can tunnel between the two wells. For N bosons the underlying Hilbert space is 
spanned by the [N + l)-dimensional Fock base {|0, N), |1, iV — 1), |7V, 0)}, where / bosons are in one 
well and iV — Z in the other well [l3l - [l^ . The initial state is prepared as a Fock state, where all the 
bosons are in one of the two wells (i.e. |0, N) or |iV, 0)), while the tunneling between the wells is turned 
off. To start the evolution, a "quench" is provided by switching on the tunneling between the two wells. 
This is realized by a sudden reduction the potential barrier between the wells in an atomic system [l[ 
or by connecting the two microwave cavities with an optical fiber @, [l^ [l7| . A similar experiment was 
performed with two atomic clouds, subject to weak interaction and separated by an adjustable potential 
barrier [ll,[l9|. 

On the theoretical side, mean-field descriptions of the Bose-Hubbard model, such as a Hartree ap- 
proximation or the Gross-Pitaevskii equation, may work well for clouds with many bosons and weak 
boson-boson interaction [T^ . However, they provide a rather poor approximation for the dynamics of 
small many-body systems (cf. Ref. [IjI)- This was also observed in a recent study by Streltsov et 
al. who compared the results of a simple Hartree (Gross-Pitaevskii) approximation with a sophisti- 
cated multi-orbital Hartree approximation 20] . The latter reveals that the bosonic clouds are related to 
superpositions of Fock states in the form of NOON states 

jiVOOiV) ^ ^ [|G, N) + e^-^^liV, 0)] . (1) 
v2 

In the following we will study the Hubbard dynamics of bosons in a double well in more detail. In 
particular, we are interested in the connection of spectral properties and the formation of NOON states, 
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based on a Fock state with all the particles in one well as the initial state. To avoid problems with 
uncontrolled approximations, we will rely on a full quantum calculation. An exact solution is available 
in a Fock-state base, as described previously in Refs. (ft, 8]. 

The paper is organized as follows: In Sect. |ll]the model, based on the Bose-Hubbard Hamiltonian, is 
defined and in Sect. Ill Al the dynamics of an isolated quantum system is explained. Then we discuss the 
dynamics of a noninteracting Bose gas in Sect. IIIII and the dynamics of an interacting Bose gas in Sect. 
IIVI The latter is divided into a study of a symmetric double well (Sect. IIV Ap and of an asymmetric 
double well (Sect. IIV Bl) . Finally, we summarize the results of our calculation in Sect. IVl and discuss 
them in Sect. EH 



II. MODEL 

The many-body Hamiltonian H of N bosons with mass m reads 

N 

+ ^t/(r„rfe), (2) 

3,k=l 

where pj is the momentum of a boson, V{rj) is the one-body potential of the double well and U is the 
two-body interaction potential. For the latter we assume that it decays very quickly with the distance 
\rj — Tfcl of the particles. This implies that particles located in different wells do not interact with each 
other. Then the many-body Hamiltonian is expressed in Fock-state representation as 

j{N-k, fc|ri, ra, rA,)(ri, ra, ynWK^t'^, r^)(|r;, r^, r'^|iV - fc', k')d?ri ■ ■ ■ Sr'^ 



N 



7 = 1 



PI 

2m 



= {N - k,k\H\N - k',k') . (3) 

For the new Hamiltonian H, which acts in the Hilbert space spanned by the Fock base, we can use the 
Bose-Hubbard Hamiltonian with local interaction in each well as a reasonable approximation 

H = J{a\a2 + 4ai) + Ui{a\ai)^ + [72(4^2)^ , (4) 

where (^i) creation (annihilation) operators for bosons in the Fock states. H, which describes 
tunneling between the two wells and the local interaction inside the well with interaction strength Uj , gives 
us a complete quantum description of the different Fock states and their superpositions. In particular, 
we can use it to study the evolution of a Fock state to a NOON state of Eq. ^ . 



A. Evolution of isolated systems 



We consider a system which is isolated from the environment. Furthermore, we assume that the system 
lives in an (A^ -1-1) -dimensional Hilbert space. With the initial state |^'o) we can get for the time evolution 
of the state 

|*t)=e-^^*|*o) (5) 
or the evolution of the return probability |(^'o|5'f)P with the amplitude 

i^ol'ft) = {^o\e-'"'\^o) ■ (6) 
In general, the amplitude (^ij^'f) can be expressed via an integral transformation of the resolvent as 

{^,\^t) = (vfi|e-'^*|vI/o) = - H)-^\^n)e-"'dz , (7) 
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where the contour F encloses aU the eigenvalues Ej (j = 0,1,..., A^) of H. With the corresponding 
eigenstates \Ej) the spectral representation of the resolvent is a rational function: 

*i(,-Hn*.).i:<iiEi<^.^, (8) 

~^ Qn+i(z) 
where Pn{z), Qn+i{z) are polynomials in z of order N , respectively, with the common denominator 

JV 



Qn+i{z) = X{{z-E,) 



These polynomials are readily evaluated by the recursive projection method (RPM) 

The expression in Eq. ^ for l^fi) = |5'o) can be interpreted as the bosonic spectral density Pi{E) 
with respect to the state |^'o): 

p,{E) i/m(*ol(i? - - - i^)-^l*o> = - E . (9) 

The amplitude of the return probability then reads as the Fourier transform of the spectral density 



(^'o|*f) = lim / p,{E)e-''^*dE . (10) 
Analogously, the overlap reads in terms of the resolvent 

{^i\-^t) = i lim / hn{^i\(E -ie- Hy^\^o)e~'^*dE (11) 

TT e^O J 



with 



\imIm{^,\{E-ie-H)-'\^'o)^7:J2{^i\E,){E,\^o)S{E-E,) , (12) 



provided that the matrix elements are symmetric. The latter is the case for the Hubbard Hamiltonian. 

The purpose of the subsequent calculation is to determine the evolution of the Fock state under the 
influence of the Bose-Hubbard Hamiltonian of Eq. In general, this is expressed in the Fock base as 

N 

j=o 

with coefficients Cj{t) = {N — j,j\^t)- For the NOON state we only need to focus on the coefficients co{t) 
and CN{t)- 

Comparing the result in Eq. ([T3|) with the expressions in Eqs. ([ro]) . (fTT|) . (fT2|) . it turns out that the 
Fourier transform of co(t) and CN(t) arc just the imaginary parts of the matrix elements of the resolvent 



Co 



{E) - lim Im{N,0\{E ~ ie ~ H)-^\N,0) {N,0\Ej){Ej\N,0)5{E - Ej) (14) 



and 



cn{E) = - lim lm{0, N\{E - ie - H)-^\N,0) {0, N\Ej){Ej\N,0)6{E - Ej) . (15) 

j 

These two expressions will be called spectral coefficients, where co{E) measures the relative weight 
\{N,0\Ej)\'^ . Integration over the energy E gives 1 for this coefficient. The coefficient cn{E) mea- 
sures the correlation between |A^, 0) and |0,7V) due to the product {0,N\Ej){Ej\N,0). The latter is real 
for a symmetric Hamiltonian. Integration over the energy E gives for this coefficient. 
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III. DOUBLE WELL: NONINTERACTING BOSE GAS 



The Bose-Hubbard Hamiltonian has two simple hmits: The local limit J = and the noninteracting 
limit Ui = U2 = 0. In the local limit for a symmetric double well with Ui = U2 pairs Fock states |iV— fc, k), 
\k,N — k) are doubly degenerate eigenstates with energy Ek = U[{N — fc)^ + fc^]. A perturbation by 
a small tunneling term will break the degeneracy. This effect is stronger at lower energies because the 
parabolic spectrum is denser there. This agrees with a numerical study [ll]- The fact that the states 
I A^, 0) and |0, N) are very close in energy may support the formation of a NOON state. 

In the absence of particle-particle interaction the Bose-Hubbard Hamiltonian Ht (i.e. the Hamiltonian 
in Eq. ^ with Ui = U2 = 0) describes only tunneling. A straightforward calculation shows that the 
eigenstate — fc; fc) of Ht with Ht\N — k;k) — J{N — 2fc)|A^ — fc; fc) has an overlap with the Fock states 
|A,0) and \0,N) as 

{N, 0\N ~k;k)= 2-^/2 (0, N\N -k;k) = (-1)'2-^/^^/^ . (16) 
This implies that the spectral coefficients of Eqs. ([T^ and p3|) have a binomial form 

co(£;)-2-^^r)<5(ii; + J(2fc-A^)) (17) 



fc=0 



cn{E) = 5] J (-l)'=<5(i? + J(2fc - N)) . (18) 

A Fourier transformation reveals a periodic behavior of the evolutionary coefficients as 

co(t) = (A,0|e-*^*|A^,0) =cos^(Jt), CN{t) = {0, N\e~'"^N,0) = {-i)^sm^{Jt) . (19) 

Thus the evolution of the Fock state leads to a NOON state with a probability that decays exponentially 
with N. This is a consequence of the fact that for an increasing N the particles disappear in the (A -I- 1)- 
dimensional Hilbert space because there is no constraint due to interaction. 



IV. DOUBLE WELL: INTERACTING BOSE GAS 



The double well with the two Fock states \N, 0), |0, N) as possible initial states can be treated within 
the RPM. This method is based on a systematic expansion of the resolvent (^'iK^ — i?)^^|4'o), starting 
from the initial base {| A, 0), |0, A^)}. The method can also be understood as a directed random walk in 
Hilbert space. This means that in comparison with the conventional random walk the directed random 
walk of the RPM visits a subspace 'H2j only once and never returns to it. In terms of N bosons, distributed 
over the double well, the subspace H2j is spanned by the base {\N — \j, N — j)}. A step from H2j to 
'H2j+2 is given by the Hamiltonian H in such a way that 'H2j+2 is created by acting H on 'H2j (cf. App. 
A). This step is provided by the tunneling of a single boson. Thus, the directed random walk follows a 
path with increasing numbers j. The directed random walk is the main advantage of the RPM which 
allows us to calculate the matrix elements (^'oK^ — -f^)~^|^o), ('^'iK-^ ~ ^^)^^|^o) of the resolvent on a 
{N + l)-dimensional Hilbert space exactly. 



A. Symmetric double well 



Now we choose Ui = U2 = U for the Bose-Hubbard Hamiltonian. Assuming that N is even, all 
projected spaces 'H2j are two-dimensional and spanned by {|A — \j,N — j)} {j ~ 0, A/2). This 
leads to a recurrence relation in the base of the two Fock states (|A, 0), |0, A)) as initial states. The 
recurrence relation reads (App. A) 
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with coefficients 



z - fk+i - J^ak{N/2 + k + l)(iV/2 - k) 



z - fk+i - .Pak{N/2 + k + l)(N/2 - k) - J^bl{N/2 + k + l)^{N/2 - k) 



(21) 



h+i 



J^bk{N/2 + k + l){N/2-k) 



and 



z - Jk+i - .PakiN/2 + k + l){N/2 - k) - ,Pbl{N/2 + k + l)^{N/2 - k) 



fk+i = UiN/2 + fc + 1)2 + U{N/2 -k-lf 



The iteration terminates after N/2 steps with 



where 



and 



9n/2 = [ h „ 

* On/2 aN/2 



aN/2 = {N,0\{z~H)-^\N,0) = {0, N\{z - H)-^\0, N) , 



bN/2 = (0, N\{z - Hy'\N, 0) = {N, 0\{z - H)-^\0, N) 



(22) 



(23) 



(24) 



(25) 



There exists an invariance of the recurrence relation under the following simultaneous sign changes in 
Eqs. IHI) and 



z — > — z, U — >■ —U, aj — —aj, bj — ^ —bj 



(26) 



This implies that a change from a repulsive to an attractive Hubbard interaction results in a mirror image 
with respect to energy of the spectral coefficients 



coiE, U) = £o{-E, ~U), cn{E, U) = ~cn{-E, -U) 



(27) 



B. Asymmetric double well 



In the case Ui — ~U2 = U we have one more variable, namely ak, bk and Ck with the following 
recurrence relations (App. A) 



9k+l - 

with matrix elements (n — N/2) 



z + Wn{k + 1) - J2(n + k+ l){n - k)bk 



D 



k+l 



(28) 



(29) 



z - 4Un{k + 1) - J2(n + k + l){n - k)ak 



D 



k+l 



(30) 



J2(n + fc + l)(n - k)ck 



(31) 
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and with 

Dk+i [z-4Un{k + l) - J^ak{n + k + l){n - k)][z + 4Un{k + 1) - J^bkin + k + l){n - k)] 

-J'^cl{n + k + lf{n-k)^ . 

The final resuh of the iteration is 

9N/2=(l''/' l^^') , (32) 

with 

aN/2 = {N,0\iz-H)-'\N,0), br,/2 ^ {0, N\iz - H)-'\0, N) , (33) 
cn/2 ^ {N,0\iz - H)-^\0,N) ^ {0,N\{z ~ H)-^\N,0) . (34) 
V. RESULTS 

The iteration of Eqs. (I^Tj) . ([^^ for a symmetric double well and the iteration of Eqs. (P^ - (PT|) for 
an asymmetric double well gives us, according to Eqs. (f24| . (|25l) and (fSSj) . (p4| . the following matrix 
elements of the resolvent 

(TV, 0|(z - Hy'\N, 0), (0, N\iz - H)~'\N, 0) = (TV, 0|(z - i/)-i|0, N) . 

These are rational functions of z, as shown in Eq. ([T])- For N bosons these are lengthy expressions with 
N + 1 poles. Therefore, it is convenient to present the results as plots with respect to energy. Examples 
of the spectral coefficients co(£') and cn{E) are shown for a symmetric double well with 100 bosons in 
Fig. [1] and with 20 bosons in Fig. [2l and for an asymmetric double well with 100 bosons in Fig. |3l A 
larger number of bosons shows a richer spectral structure. The diagonal coefficient co{E) in the case of 
100 bosons is remarkably different from the off-diagonal coefficient cn{E) because the latter does not 
have spectral weight from eigenstates whose energy Ej is larger than the energy of the initial Fock state 
E = UN'^. The reason for this feature is the double degeneracy of the eigenvalues mentioned in Sect. HID 
The signs of the product (0, N\Ej){Ej\N, 0) for adjacent eigenvalues are opposite to each other. Since the 
eigenvalues get closer pairwise as we increase their energy, the contribution of the two levels cancel each 
other for each pair inside the sum of Eq. ([15]). This interaction effect is also visible for 20 bosons (Fig. 
[5]), although the cancellation is incomplete then due to a larger level distance. This can be considered as 
an effect of spectral fragmentation, where the spectrum has a nondegenerate low-energy and a degenerate 
high-energy part, caused by the competition of tunneling and interaction. 

The contribution of the two Fock states |0, A^), \N,0) to the evolution in Eq. (fT3)) is given by the 
coefficients co(i) = (iV, 0|4't), CAr(i) = (0,iV|4't). In Fig. |3]the real parts of these coefficients are plotted 
for 100 bosons. Their evolution indicates a collapse and revival behavior. The plot of |co(i)|, |cjv(OI 
as a two-dimensional vector for 20 bosons in Fig. [S] shows a complex dynamical behavior that cannot 
be described by a simple equation of motion. This observation suggests a statistical description with a 
probability P(|co(i)|, |cAr(t)|) which measures how often certain values of |co(i)|, |cAr(i)| are visited during 
the evolution in a period of time. The result for 20 bosons is plotted in Fig. [6]for [/ = 0.1, J = VT5. 
It indicates that there is a strong correlation between the coefficients, where the most favored values are 
\co{t)\ « \cN{t)\- 0.35. 

For an asymmetric double well with interaction strength ±U the spectrum is different because of 
the absence of double degeneracy of the eigenvalues (cf. Fig. [3]). There are two "bands", one around 
E = UN'^, the other around E = —UN'^, where the widths of the bands is characterized by the tunneling 
rate J. Moreover, the off-diagonal part cn{E) appears closer to zero energy and its values are very small. 
This indicates that the off-diagonal part has overlaps with energy levels which are different from those 
of the diagonal part co{E) . For the evolution only the latter contribute substantially, preventing the 
system to create a NOON state. 
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VI. DISCUSSION AND CONCLUSIONS 

In order to understand the evolution of an isolated many-body bosonic system, we start with noninter- 
acting bosons (i.e. Ui =U2 = 0) of Sect. IIIII The spectral properties are characterized by (i) equidistant 
energy levels with distance J and (ii) a binomial weight distribution of the energy levels. The evolution 
of a Fock state is characterized by a periodic behavior with a single frequency w = J/27r as a direct 
consequence of the equidistant energy levels. The amplitudes for visiting the initial Fock state |A^, 0) or 
the complimentary Fock state |0, N) vary with cos^(Jt) or {—i)^ sin^(Jt), respectively. This implies for 
a large number TV of bosons that (i) these states are visited only for a very short period of time and (ii) 
the two Fock states are visited at different times. Thus the formation of a NOON state is very unlikely 
for noninteracting bosons. 

A simple qualitative picture for the general evolution of the Fock state is the random walk in Hilbert 
space. In case of noninteracting bosons the particles can walk independently of each other which enables 
them to explore the entire Hilbert space spanned by the Fock states without restriction. A simultaneous 
overlap of with both Fock states |iV, 0) and |0,A^) is very unlikely then. Once we have turned on 
the boson-boson interaction the particles experience a mutual influence which restricts their individual 
random walks. This is related to the fact that the systems stays much longer in the energetically (almost) 
degenerate Fock states \N,0) and |0,A^) than in the noninteracting case (cf. Fig. [S]) and, what is even 
more important here, they can have a simultaneous overlap with both Fock states, such that they create 
a NOON state. In terms of the spectral properties the interaction modifies (i) the energy levels, which 
are not equally spaced, and (ii) the weight distribution of the levels, which are not binomial any longer 
(cf. Fig. [T][3]). This, of course, affects also the evolution of the Fock state which is more complex now, 
since many different frequencies are involved. A particular feature is the spectral fragmentation (cf. Fig. 
[T]), where only a part of the spectrum contributes to the off-diagonal coefficient cn{E). This is a kind of 
Hilbert-space localization, where transitions to the high-energy part of the Hilbert space are completely 
suppressed, similar to the self-trapping found in the Hartree approximation of the Bose-Hubbard model 
[l3|. It should be noticed, however, that spectral fragmentation appears at a much weaker interaction 
strength than the self-trapping effect. For t/ ~ J, which is the threshold for self-trapping [l^, there is 
only one eigenvalue with significant weight |(7V, 0|£^j)p @. Thus it is not clear whether or not the two 
effects are directly are connected. 

For the asymmetric double well the situation is different due to the existence of two "bands" and the 
absence of the double degeneracy. The main consequence is the absence of a support for the formation of 
NOON states because the off-diagonal coefficient CAr(i) is strongly suppressed. From this observation we 
can conclude that the evolutionary entanglement is much more favorable in the symmetric double well. 
This is in agreement with the results of the multi-orbital Hartree calculation of Ref. 20] . 

In conclusion, we have studied the evolution of a bosonic Fock state |A^, 0) in a double well and 
found that a local particle-particle interaction supports the formation of a NOON state, provided that the 
interaction is not too strong. This is accompanied by a fragmentation of the spectrum. The latter is 
characterized by the fact that only eigenstates with energies less than the energy of the initial Fock state 
can be reached in the evolution. This interaction effect causes a Hilbert-space localization and prevents 
the evolution of the Fock state to disappear in the depth of the Hilbert space. This is the main reason for 
a favorable creation of a NOON state. The appearance of a NOON state is suppressed though for strong 
interaction because then the restriction of the Hilbert space is too severe and does not allow to reach the 
complementary Fock state |0, A^). 
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Appendix A: Recursive projection method 

Given is a sequence of projectors Pj {n > j > 0), defined by the recurrence relation 

P2k+i = P2k-i -P2k {n>k> 0) 
with initial conditions P_i = 1, Pq said by the Hamiltonian H through the properties 

P2kHP2k+l = P2kHP2k+2, P2k+lHP2k - P2k+2HP2k ■ (Al) 

The projection of the resolvent [z — H)^^ defines 

5fc = P2{n-k) {Z - ff)2(!i-/c)-l P2{n~k) (0 < fc < n) , (A2) 

where (•■•)2(ri-fc) inverse on the P2(n-fe)~projected Hilbert space. Then satisfies the recurrence 

relation 

9k = {z- hk)2ln-k) (A3) 



, _ , P2nHP2n k^O 

" ^ P2{n-k)HP2(n-k) + P2{n-k)H gk-lH P2(n-k) l<k<n ' ^"^^^ 



with 



Of interest is here only the case k — n, where we have from Eq. (IA2p 

gn^Po{z~H)-'Po . 
For the specific case of the double well we choose n = N/2 and the projectors 

Po = |7V,0)(7V,0| + |0,7V)(0,iV|, P2 = |7V-l,l)(7V-l,l| + |l,iV-l)(l,iV-l|,..., 

Pat \N/2,N/2){N/2,N/2\ . 

With the Hubbard Hamiltonian of Eq. ^ the diagonal terms of the effective Hamiltonian in Eq. (jA4 
read 

P2(n-k)HP2(n-k) ^ Pn -2kP[ Pn ~2k 



= [Ui{Nl2 + kf + U2{N/2 - \N/2 + k, N 12 - k) {N/2 + k, N/2 - k\ 
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FIG. 1: Spectral coefficients of Eqs. (O and ^ for 100 bosons with U/J « 0.023 and e = 0.01. The energy 
of the initial Fock state is -B = 1000. The spectral fragmentation appears around E « 970, where the levels 
are nondegenerate at lower energies but almost degenerate for higher energies. This is a consequence of the 
competition between tunneling and interaction, which the latter wins at higher energies. 



+ [Ui{N/2 - kf + U2{N/2 + kf] \N/2 - fc, N/2 + k) {N/2 - k, N/2 + k\ 
The off-diagonal terms of the effective Hamiltonian in Eq. (|A4[) read 



P2{n-k)HP2(n-k+l) — P2N -2kH Pn -2k+2 



-JyjN/2 + k^jN/2 - /c + 1 



{\N/2 + k, N/2 - k){N/2 + - 1, N/2 - fc + 1| + \N/2 - k, N/2 + k){N/2 -k + I, N/2 + fc - 1|) 
This leads for Ui = f72 to Eqs. ([21]), ([22]) and for J7i = -J/z to Eqs. ([291), jSQl) and 
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FIG. 2: Spectral coefficients for 20 bosons with U/J~ 0.1 and e = 0.01. The energy of the initial Fock state is 
E = 40. The almost degenerate states appear above 42. 



11 




1 .5e-040 



1e-040 



5e-041 



-5e-041 



-1e-040 



1500 



1500 



FIG. 3: Spectral coefficients for 100 bosons in an asymmetric double well with Uj J « ±0.023 and e = 0.01. The 
energy of the Fock states |0, A^) and \N,Q) is E = ±1000, respectively. 
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FIG. 4: Evolution of the real part of the evolutionary coefficients co(t) = (A'^, 0|'I't) (upper panel), CAr(t) = 
(0, Nl'i/t} (lower panel) for 100 bosons with U / J ^ 0.023. The time scale is given in inverse units of the interaction 
strength O.lh/U and the energy of the initial state is .E = 1000. 
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trajectory (|co(t)|, |cjv(t)|) starts at (1,0). 




FIG. 6: Distribution P{\co{t)\, |cjv(t)|) of \co{t)\, cjv(t)j over a time period of OAh/U for 20 bosons with U/J ^ 
0.026. The axes are scaled by a factor 100 and the vertical axis is in arbitrary units. This plot indicates a strong 
correlation between the two spectral coefficients, supporting the formation of a NOON state. 



